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ABSTRACT 


A  parametric  study  of  the  speeds  of  the  current  sheet  and  center-of- 
gravity  of  current  in  a  magnetic  annular  shock  tube  has  been  carried  out. 
The  parameters  varied  include:  (1)  the  drive  current,  (2)  the  polarity  of 
the  inner  electrode,  (3)  the  gas  pressure,  (4)  the  nature  of  the  gas,  (5)  the 
radius  of  the  inner  electrode  and  (6)  the  material  of  the  inner  electrode. 

The  most  interesting  result  of  the  investigation  was  the  observation  of 
limiting  speeds  for  the  current  sheet  and  center-of-gravity  of  current  of 
approximately  8  cm/jU.s  and  3  cm/jLi.s,  respectively.  These  speeds  were 
not  exceeded  even  under  conditions  where  the  magnetic  pressure  exceeded 
the  dynamic  pressure  by  a  factor  of  10.  The  most  probable  explanation  of 
the  limiting  speed  is  that  it  is  due  to  the  inertial  drag  of  material  ablated 
from  the  insulator  at  the  driver  end  of  the  shock  tube. 


INTRODUCTION 


In  a  magnetic  annular  accelerator,  gas  in  the  space  between  two 
coaxial  electrodes  is  accelerated  by  the  magnetic  pressure  behind  a  cur¬ 
rent  sheet  maintained  between  the  electrodes  by  a  condenser  discharge. 

Such  accelerators  have  potential  applications  as  engines  for  space  propul¬ 
sion,  plasma  injectors  for  fusion  devices,  and  shock  tubes  for  the  study 
of  gas  dynamics  in  the  speed  range  beyond  that  which  can  be  reached  by 
thermal  or  chemical  means.  Where  the  interest  is  in  shock  tube  appli¬ 
cations,  the  device  has  been  called  a  "magnetic  annular  shock  tube"^  and 
the  gas  is  introduced  at  a  uniform  pressure;  where  the  interest  is  in  pro- 
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pulsion  or  plasma  injection,  the  accelerator  is  often  called  a  "plasma  gun" 
and  the  gas  is  introduced  as  a  puff.  Although  considerable  study  has  been 
devoted  to  both  modes  of  operation,  neither  practical  applications  nor  a 
satisfactory  theory  of  the  operation  of  such  accelerators  have  yet  been 
achieved. 

The  present  investigation  represents  an  attempt  to  obtain  a  better 

understanding  of  the  factors  effecting  the  speed  of  the  current  sheet  in  a 

magnetic  annular  shock  tube  having  a  large  radius  ratio  and  no  externally 

applied  bias  fields.  It  was  motivated  by  evidence  obtained  in  a  previous 
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investigation  of  the  shape  of  the  current  sheet  that  serious  discrepancies 
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between  the  measured  and  predicted  ’  sheet  speeds  were  occurring.  Since 
the  sheet  speed  is  one  of  the  most  important  and,  in  principle,  easiest  param¬ 
eters  to  predict,  this  was  considered  quite  serious.  In  this  connection  it 
should  be  noted  that  in  a  similar  shock  tube  having  a  small  radius  ratio  and 
magnetic  bias  fields,  Patrick^  has  produced  magnetically  driven  shock  waves 
which  do  move  at  the  theoretical  velocity.  Although  the  origin  of  this  appar- 

A  Y 

ent  difference  in  operation  is  not  fully  understood,  Fishman  and  Petschek  ’ 
have  shown  that  both  the  radius  ratio  and  the  bias  fields  have  an  important 
effect  on  the  performance  of  a  magnetic  annular  shock  tube  and  that  Patrick's 


shock  tube  may  be  expected  to  operate  in  a  completely  different  manner 
than  the  one  used  in  the  present  investigation.  Thus  to  a  large  extent, 
Patrick's  results  must  be  regarded  as  applying  to  a  different  class  of  shock 
tubes,  and  there  is  no  real  contradiction  between  his  results  and  those  re¬ 
ported  in  this  paper. 

The  usual  approximation  made  in  obtaining  an  expression  for  the 
sheet  speed  is  to  use  the  "slug  model.  "  This  model  assumes  that  all  the 
gas  swept  up  by  the  current  sheet  moves  at  the  sheet  speed  and  that  wall 
friction  may  be  neglected.  Under  these  conditions  the  momentum  equation 
becomes  (MKS  units) 


Pi  dz  + 


P^U 


(l-p^/p^) 


(1) 


where  B  is  the  magnetic  field  immediately  behind  the  current  sheet,  is 
the  vacuum  permeability,  u  is  the  speed  of  the  current  sheet,  z  is  the  axial 
position  of  the  front  of  the  slug,  and  Pj^  and  are,  respectively,  the  gas 
densities  ahead  of  and  in  the  slug.  In  order  to  obtain  the  density  ratio  P^/P2 
one  must  employ  the  equations  of  mass  and  energy  conservation  and  the 
equation  of  state.  However,  in  most  of  the  cases  we  shall  consider  P1/P2 
small  enough  so  that  it  can  be  neglected  compared  to  one  without  significant 
error . 

Equation  (1)  provides  a  fundamental  check  on  the  operation  of  any 

magnetic  accelerator.  Failure  of  the  measurements  to  fulfill  this  equation 

implies  a  basic  difficulty  in  either  the  operation  of  the  accelerator  or  the 

experimental  techniques.  Unfortunately,  it  is  difficult  to  use  this  equation 
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in  the  case  of  plasma  guns,  ’  ’  where  the  gas  is  introduced  in  a  puff,  due 

to  the  uncertainty  in  the  initial  density  distribution.  Thus,  if  achieving  an 
understanding  of  the  acceleration  mechanism  is  the  goal  of  an  experiment, 
operation  of  the  accelerator  as  a  shock  tube  at  constant  initial  density  is 
desirable.  A  number  of  such  shock  tube  investigations^’  have 

been  carried  out  and  although  some  measure  of  agreement  with  Eq.  (1)  has 
been  reported,  the  situation  is  far  from  satisfactory.  Most  of  the  measure¬ 
ments  have  been  made  in  shock  tubes  having  a  radius  ratio  greater  than  3  so 


-2- 


that  the  magnetic  pressure  varies  by  more  than  a  factor  of  9  across  the 

annulus  and  there  exists  considerable  ambiguity  about  just  how  one  should 
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employ  Eq.  (1).  Furthermore,  it  has  been  shown  ’  '  that  the  sheet  speed 
in  a  shock  tube  of  typical  design  depends  on  the  polarity  of  the  inner  electrode. 
Finally,  in  most  cases  the  current  pulse  had  the  form  of  a  damped  sine  wave 
so  that  the  speed  varied  continuously  and  was,  therefore,  difficult  to  deter¬ 
mine  precisely.  These  considerations  combined  with  the  fact  that  the  results 
of  different  investigators  employing  similar  devices  do  not  agree  with  the 
precision  one  would  expect,  have  produced  a  confused  situation  which  leaves 
one  in  considerable  doubt  as  to  what  is  occurring  in  such  large  radius  ratio 
magnetic  annular  shock  tubes. 

In  the  present  study  we  have  concentrated  on  measurements  of  the 
speed  of  the  current  sheet  and  back  EMF  in  the  shock  tube.  To  provide  in¬ 
formation  which  would  give  a  reasonable  check  on  the  predictions  of  the 
momentum  equation  and  provide  an  overlap  with  the  work  of  other  investi¬ 
gators,  as  many  of  the  parameters  effecting  operation  of  the  shock  tube 
were  varied  as  possible  within  the  physical  limitations  of  the  experimental 
setup.  These  include:  (1)  the  driving  current,  (2)  the  polarity  of  the  inner 
electrode,  (3)  the  gas  pressure,  (4)  the  nature  of  the  gas,  (5)  the  radius  of 
the  inner  electrode,  and  (6)  material  of  the  inner  electrode.  The  operation 
was  also  studied  with  the  center  electrode  replaced  by  a  button  to  simulate 
conditions  encountered  in  "conical"  shock  tubes.  The  only  parameters  not 
varied  were  the  radius  of  the  outer  electrode  and  the  material  of  the  insula¬ 
tors  at  the  ends  of  the  device.  As  we  shall  see  in  the  discussion,  failure 
to  vary  the  material  of  the  insulators  deprived  us  of  a  valuable  piece  of 
information  bearing  on  the  interpretation  of  the  results,  but  the  shock  tube 
was  dismantled  before  this  could  be  easily  rectified. 

The  most  significant  result  of  the  investigation  was  the  observation 
of  a  "limiting  speed"  beyond  which  it  was  not  possible  to  drive  the  current 
sheet  even  under  conditions  where  the  total  magnetic  pressure  exceeded 
the  dynamic  pressure  by  more  than  an  order  of  magnitude.  This  "limiting 
speed"  was  substantially  independent  of  any  of  the  parameters  which  were 
varied  in  the  experiment  and  provides  clear  evidence  for  an  important  source 
of  drag  on  the  current  sheet  in  addition  to  that  due  to  the  dynamic  pressure  of 


-3- 


the  gas  initially  in  the  shock  tube.  Estimates  of  the  effect  of  the  boundary 
layer  on  the  walls  of  the  shock  tube  show  that  it  is  unlikely  to  be  the  result 
of  wall  friction.  It  is  also  unlikely  to  be  associated  with  the  drag  due  to 
positive  ions  striking  the  cathode  in  the  case  where  the  ions  carry  the  cur¬ 
rent.  The  most  plausible  hypothesis  is  that  it  is  associated  with  the  inertia 
of  material  ablated  from  the  insulator  at  the  driving  end  of  the  shock  tube. 

This  material  is  ionized  either  by  photons  or  electron  impact  and  becomes 
entrained  by  the  currents.  A  simplified  theory  of  the  process  based  on  the 
assumption  of  a  steady  state  ablating  boundary  layer  on  the  insulating  wall 
leads  to  a  qualitatively  correct  description  of  the  observations.  To  make 
the  theory  quantitative,  one  must  calculate  the  thickness  of  the  boundary 
layer.  This  is  a  problem  beyond  the  scope  of  the  present  paper  which 
involves  (1)  the  heat  of  vaporization  and  vapor  pressure  law  for  the  wall 
material,  (2)  the  electrical  and  thermal  conductivity  of  the  ablation  products, 

(3)  cross-sections  for  ionization  of  ablation  products  by  photons  and  electrons, 

(4)  energy  dissipation  by  ion  slip  and  (5)  photon  production  in  the  gas  and  at 
the  electrodes.  Nevertheless  the  problem  does  appear  tractable  and  a  full 
theory  of  such  a  boundary  layer  should  have  many  important  applications. 

EXPERIMENTAL  APPARATUS 

A  schematic  diagram  of  the  magnetic  annular  shock  tube  used  in  the 

experiments  is  shown  in  Fig.  1.  It  is  identical  with  that  used  in  the  previous 
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studies  of  shock  shape  and  consists  of  a  cylindrical  copper  outer  electrode 
6"  in  diameter  and  34"  long  and  an  interchangeable  inner  electrode.  The 
insulators  at  the  ends  of  the  device  were  Incite  sealed  to  the  electrodes  by 
O  rings.  The  tube  was  pumped  through  a  series  of  8  slots  arranged  symmet¬ 
rically  around  the  outer  electrode  2"  from  the  driver  end  and  could  be  evacu¬ 
ated  to  a  pressure  of  a  few  tenths  of  a  micron. 

To  preserve  cylindrical  symmetry,  the  shock  tube  was  connected  to 
the  condenser  bank  through  a  2'  section  of  coaxial  line  terminating  in  a 
triggered  spark  gap.  The  measured  inductance  L^  of  this  line  was  .  06  fjh. 

The  condenser  bank  was  connected  to  form  a  lumped  constant  line  having  a 
total  capacity  of  384  fji  and  a  characteristic  impedance  of  .  023  ohms. 
The  system  operated  reliably  at  positive  or  negative  voltages  between 
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2  and  10  KV  and  produced  a  current  pulse  which  rose  with  a  time  constant 
of  1 .  5  jU  sec  and  was  approximately  constant  for  about  10  fj.  sec  when  dis¬ 
charged  through  a  matched  load.  Under  most  conditions  the  impedance  of 
the  shock  tube  was  small  compared  tO'  and  the  steady  value  of  the  current 
was  approximately  V^/Z^. 

To  maintain  reasonable  gas  purity  a  continuously  operating  flow  sys¬ 
tem  was  used  which  changed  the  gas  in  the  tube  once  a  minute.  The  pressure 
was  measured  by  an  Alphatron  gauge. 

Provision  for  preionizing  the  gas  was  provided  by  four  small  button 
electrodes  symmetrically  placed  in  the  insulating  wall  at  the  driver  end  of 
the  shock  tube.  These  were  individually  connected  to  four  2  jUf-10  KV  con¬ 
densers  which  could  be  discharged  simultaneously  through  a  single  switch. 
Preionization  made  no  significant  difference  in  the  operation  of  the  shock 
tube  after  breakdown,  but  it  did  eliminate  delays  in  breakdown  and  was 
normally  employed  for  this  reason.  It  may  be  noted  in  this  connection  that 
no  matter  where  the  initial  breakdown  occurred  in  the  tube,  the  main  cur¬ 
rent  sheet  always  originated  at  the  driving  end.  The  reason  for  this  is  that 
radiation  from  the  initial  streamer  photoionizes  the  gas  and  makes  it  a 
sufficiently  good  conductor  to  prevent  further  penetration  of  the  magnetic 
field. 

The  "optical  speed"  of  the  discharge  was  measured  by  a  series  of 
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5  tungsten  photoelectric  detectors  of  a  type  developed  by  Camac  spaced 

at  6"  intervals  along  the  length  of  the  tube.  Tungsten  detectors  operating  in 
UV  were  used  in  preference  to  photomultipliers  operating  in  the  visible  be¬ 
cause  they  gave  cleaner  signals,  are  extremely  stable  and  are  uneffected  by 
scattered  light.  The  light  entering  the  tungsten  detectors  was  collimated 

by  two  1  mm  holes  spaced  2  cm  apart  giving  a  solid  angle  area  product 
-5  2 

S2A=1.5xl0  cm  ster  and  a  maximum  spot  diameter  of  7.  5  mm  at  the 
axis  of  the  tube. 

The  current  wave  form  was  measured  by  a  one  turn  rectangular  coil 
feeding  directly  into  a  terminated  93  cable.  Voltage  wave  forms  were  meas¬ 
ured  through  100  K  attenuators  connected  to  terminals  at  both  ends  of  the 
center  electrode  as  indicated  in  Fig.  1.  An  additional  flux  loop  passing 
through  the  center  electrode  and  around  the  test  section  was  used  as  a  check 
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on  lon.g  3-s  the  driver  currents  do  not  penetrate  to  the  outside  of 

the  electrodes  the  voltage  from  the  flux  loop  should  be  identical  to 

This  was  always  the  case  in  the  present  experiments. 

To  study  conditions  in.  the  interior  of  the  shock  tube  scanning  probes 
could  be  inserted  through  O-ring  seals  in  the  outside  electrode  or  the  terminal 
wall  as  shown  in  Fig.  1. 

EXPERIMENTAL  MEASUREMENTS 
Summary  of  Conditions 

Measurements  were  made  over  a  range  of  driver  currents  from 
40  to  400  ka  and  initial  pressures  p..  from  50  to  1000  |i  Hg  corresponding  to 
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initial  particle  densities  n^^  from  1.8x10  '  to3.6xl0  cm  .  The  test 
gases  employed  were  hydrogen,  helium,  air  and  argon.  Copper  center 
electrodes  l/4",  1",  Z"  and  4"  in  diameter  and  a  stainless  steel  electrode 
Z"  in  diameter  were  used  in  the  study.  In  addition,  a  short  series  of  ex¬ 
periments  was  performed  in  which  all  but  a  short  segment  of  the  center 
electrode  was  omitted  or  replaced  by  teflon  or  pyrex  rods. 

In  all  cases  the  center  electrode  was  run  both  positive  and  negative 
3 

since  previous  experiments  had  indicated  that  the  shape  and  speed  of  the 
current  sheet  were  polarity  dependent.  Although  the  present  measurements 
also  showed  a  polarity  dependence  the  effects  were  relatively  small,  and  we 
have  in  general  given  the  results  only  for  the  negative  case  since  this  has 
been  the  most  common  mode  of  operation  in  other  laboratories. 

Current  and  Voltage  Measurements 

A  typical  set  of  current  and  voltage  wave  forms  for  both  positive  and 
negative  center  electrode  is  shown  in  Fig.  Z.  The  traces  from  top  to  bottom 
are  the  flux  loop  voltage  the  driving  voltage  the  driver  current  I^ 

and  the  terminal  voltage  R-sasonably  steady  conditions  are  obtained 

for  a  test  time  of  10  to  15  /i  sec  and  except  for  the  polarity  of  the  signals, 
there  is  relatively  little  difference  between  positive  and  negative  operation. 

It  is  interesting  to  note  that  while  and  ^^2  exceed  IZOO  V  during 

the  test  time,  V^z  less  than  50  V  until  the  current  interface  strikes 

the  terminal  insulator  at  which  time  it  rises  abruptly  to  a  value  of  over  300  V. 
This  shows  that  only  a  small  fraction  of  the  total  power  supplied  to  the  device 
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is  being  dissipated  by  the  arc  drop  and  that  the  gas  ahead  of  the  current 

sheet  is  a  sufficiently  good  conductor  to  prevent  any  appreciably  diffusion 

of  current  to  the  terminal  wall.  It  is  also  interesting  to  note  that  while 

V .  r  1  V,  and  I  oscillate  as  a  result  of  the  impedance  mismatch  between 
45  12  o  ^ 

the  shock  tube  and  the  current  source,  ^ ^2  changes  sign.  Since  ^2 

has  the  sign  of  the  current  at  the  terminal  end,  this  shows  that  current  loops 
are  being  formed  within  the  shock  tube  each  time  the  driver  current  changes 
sign. 

An  important  parameter  which  may  be  obtained  from  the  type  of  data 
shown  in  Fig.  2  is  the  flux  speed  u^^.  The  flux  speed  is  the  axial  speed  of 
the  center-of-gravity  of  the  current  density  in  the  shock  .tube  and  is  defined 
by  the  equation 


where  d^/dt  =  is  the  rate  of  change  of  magnetic  flux  through  the  flux 
loop,  L'  =  (jL(,^/2  w)  I  n  (r^/r^;)  is  the  inductance  per  unit  length  of  the  shock 
and  (r2/i'j^)  is  radius  ratio  for  the  electrodes.  During  the  steady  test  time 
the  term  involving  dl/dt  may  be  neglected  and 

where  E  .-  =  V.^/r  jl  n  (r-/r, )  and  B  =  U  I  /2  tt  r.  It  should  be  noted  that 
is  not  the  E/B  velocity  at  the  driver  end  except  when  the  terminal 
voltage  V22  is  identically  zero  and  that  measurement  of  the  driving  voltage 
V,  ,  alone  is  insufficient  to  determine  whether  there  is  any  axial  motion  of 
the  current  in  a  coaxial  accelerator. 

A  plot  of  u,  as  a  function  of  the  parametric  speed 

i-J 

r  2  /  I 

ui  =  [b"^  (4) 

is  shown  in  Fig.  3  for  a  variety  of  operating  conditions.  Note  that  u^  is  the 
speed  obtained  from  Eq.  (1)  by  balancing  dynamic  and  magnetic  pressures 
at  the  inside  electrode  and  neglecting  the  term  Also  shown  in  Fig.  3 

by  the  diagonal  lines  are  the  parametric  speeds 
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(5) 


■  r  2  /  1 

'"Z  =  L®  o  j 

obtained  by  balancing  dynamic  and  magnetic  pressures  at  the  outside  elec¬ 
trode,  and 

u  =  [l/2  L'l^  /,r  (r^  -  '  (6) 

obtained  by  balancing  the  integrated  dynamic  and  magnetic  pressures.  When 
plotted  as  a  function  of  all  of  the  data  fall  more  or  less  on  a  common  curve 
which  is  asymptotic  to  Uj^  for  small  values  of  u^,  but  levels  off  at  a  value’  of 
3  cm/jX  s  for  large  values  of  Uj^.  When  plotted  as  a  function  of  either  u^  or  u 
no  simple  correlation  is  evident  and  the  data  scatter  widely.  This  shows 
that  for  high  values  of  u^^  Eq.  (1)  is  not  valid  and  there  must  exist  an  impor-. 
tant  drag  on  the  current  in  the  shock  tube  in  addition  to  the  inertial  drag  of 
the  test  gas. 

A  second  parameter  which  may  be  obtained  from  the  type  of  data 

shown  in  Fig.  2  is  the  value  of  terminal  voltage  Typical  values  of 

^22  steady  test  time  are  shown  in  Fig.  4  as  a  function  of  the 

current  per  unit  circumference  of  the  inner  electrode  1^/2  Trr^^.  This  gave 

the  best  simple  correlation  of  the  data  that  could  be  found.  It  can  be  seen 

that  y 22  roughly  inversely  proportional  to  I^/ZTtrj^.  If  we  assume  that 

y 22  related  to  the  arc  drop  at  the  center  electrode  and  that  the  arc  spot 

has  a  length  of  about  10  cm  as  indicated  by  data  to  be  presented  later,  this 

leads  to  an  approximately  constant  power  dissipation  in  the  boundary  layer 

at  the  inner  electrode  of  40  Kw/cm  .  If  this  power  goes  into  the  copper  it 

will  only  raise  the  surface  temperature  of  the  wall  of  the  order  of  10°K 

which  is  insignificant.  However,  if  it  goes  into  ablation  of  a  surface  layer 

of  impurities,  as  it  probably  does,  it  is  sufficient  to  produce  a  particle  flux 
23  -1  -2 

of  the  order  of  10  sec  cm  which  is  comparable  to  the  particle  flux 
incident  on  the  current  sheet  and  could  have  a  significant  effect  on  the  flow 
near  the  electrodes. 

Ultraviolet  Light  Measurements 

A  t'j'pical  set  of  oscillograms  obtained  using  the  tungsten  U'V  photo¬ 
electric  detectors  is  shown  in  Fig.  5  for  discharges  in  H^,  H^,  Air  and  Ar 
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all  at  the  same  pressure  and  driving  voltage.  The  lower  traces  show  the 

outputs  of  the  five  UV  detectors  feeding  into  a  terminated  93  cable  and 

the  upper  traces  show  the  terminal  voltage  ^^2  which  gives  an  indication 

of  the  arrival  time  of  the  current  sheet  at  the  end  wall  of  the  shock  tube  4" 

beyond  the  last  detector.  By  inserting  probes  through  the  terminal  wall 

of  the  shock  tube  to  shadow  the  UV  detectors,  it  has  been  established  that 

most  of  the  light  observed  comes  from  the  arc  spot  on  the  inside  electrode. 

It- has  further  been  established  that  the  position  of  this  spot  correlates 

exactly  with  the  position  of  the  leading  current  sheet  in  the  shock  tube  and 

.3 

that  no  significant  shock  wave  precedes- it.  Thus,  the  records  in  Fig.  5 
give  the  speed  of  the  leading  current  sheet  down  the  shock  tube.  It  can  be 
seen  that,  except  for  the  first  interval  during  which  the  current  is  rising, 
the  sheet  speed  is  essentially  constant  over  the  entire  length  of  the  shock 
tube.  It  can  also  be  seen  that  this  speed  is  very  nearly  independent  of  the 
molecular  weight  of  the  gas,  which  again  indicates  that  the  major  drag  on 
the  current  is  not  due  to  the  inertia  of  .the  test  gas.  According  to  Eq.  (1), 
a  current  sheet  propagating  in  H2  at  uniform  speed  should  move  at  4.  5' time 
the  speed  of  a  sheet  propagating  in  Ar  when  driven  by  the  same  magnetic 
pressure.  This  is  clearly  not  the  case  and  the  observed  speed  in  is  only 
about  20%  higher  than  that  in  Ar. 

A  somewhat, more  quantitative  picture  of  the  situation  is  given  in 
Fig.  6  which  shows  the  measured  sheet  speeds  u^  plotted  as  a  function  of 
time  for  both  positive  and  negative  center  electrodes.  In  addition  to  the 
features  just  discussed,  the  data  in  Fig.  6  confirms  previously  reported 
results  that  there  is  a  polarity  asymmetry  with  respect  to  u^  in  which  the 
sheet  tends  to  propagate  slightly  faster  when  the  center  electrode  is  positive. 

Also  shown  in  Fig.  6  is  the  flux  speed  u^^  deduced  from  the  flux  loop 
voltage  Within  the  accuracy  of  measurement  this  was  the  same  in  all 

cases  and  only  a  single  curve  is  shown.  Note  that  u^^  was  constant  over  the 
entire  range  for  which  it  could  be  determined  even  though  and  V^g  were 
constant  only  in  the  interval  from  5  to  15  jj,  sec.  It  can  be  seen  that  u^^  is 
less  than  half  as  large  as.  u^  which  indicates  that  the  center-of-gravity  of 
the  current  density  is  lagging  considerably  behind  the  position  of  the  leading 
current  sheet.  This  suggests  that  estimates  of  the  specific  impulse  of  a 
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coaxial  accelerator  based  on  measurements  of  may  be  seriously  in  error, 

since  it  is  more  reasonable  to  associate  the  speed  of  the  bulk  of  the  material 

with  Ut  than  u  . 

L  s 

A.plot  of  the  sheet  speed  u^  as  a  function  of  the  parametric  speed  Uj^ 
is  shown  in  Fig.  7  for  conditions  corresponding  to  those  of  Fig.  3.  It  can 
be  seen  that  u^  exhibits  the  same  general  trend  as  the  flux  speed  u^^  except 
that  the  correlation  with  Uj^  is  not  quite  as  good  and  the  limiting  velocity  is 
not  quite  so  well  defined  since  u^  is  still  rising  slowly  at  the  highest  values 
of  Uj^  investigated. 

In  addition  to  the  speed  of  the  current  sheet,  we  may  also  obtain  an 
estimate  of  the  intensity  of  UV  radiation  in  the  shock  tube  from  photoelectric 
records  of  the  type  shown  in  Fig.  5.  The  photoelectric  current  1^  in  the 
tungsten  detectors  is  related  to  the  spectral  intensity  d(^/dnAdv  of  the  inci¬ 
dent  UV  radiation  by  the  expression: 

I  =eflA/e  (d^/dndAdv)dv  (7) 

p  V 

where  e  is  the  electronic  charge,  OA  is  the  solid  angle  area  product  for 

the  collimating  holes  and  is  the  quantum  efficiency  of  the  tungsten  cathode. 

Data  on  the  photoelectric  efficiency  of  untreated  tungsten  summarized  by 
15 

Weissler  shows  that  rises  sharply  from  a  threshold  at  hv  w  10  ev  to 
an  approximately  constant  value  of  0.  1.  Thus  to  a  first  approximation  we 
may  write 

Ip  =  0.  1  eilA  {d^/dadA)  (8) 

where  d(^i/dJ2dA  is  the  flux  of  photons  per  unit  area  per  unit  solid  angle  with 
energies  above  10  ev  falling  on  the  detector.  Since  it  is  known  that  most 
of  the  photons  entering  the  detector  come  from  the  arc  spot  on  the  inside 
electrode,  we  may  obtain  the  total  intensity  per  unit  solid  angle  d^/dSl  normal 
to  the  electrode  by  integrating  the  output  of  the  photo  detectors  with  respect 
to  time  and  multiply  the  result  by  ZTrr^u^.  This  leads  to  the  expression: 

d(/»/dn  =  20  Trrj^u^  (/  Ipdt)/e  Q  A.  (9) 
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A  plot  of  d(^/dSl  as  a  function  of  the  driving  current  is  shown  in  Fig.  8. 

As  can  be  seen  in  Fig.  5,  signals  from  the  UV  detectors  fluctuated  con¬ 
siderably  and  the  reproducibility  of  the  measurements  was  about  a  factor 
of  3.  Neveretheless,  the  order  of  magnitude  of  d(^/df2  is  certainly  correct 
and  shows  that  if  the  angular  distribution  is  isotropic,  there  is  roughly  one 
photon  with  energy  in  excess  of  10  ev  emitted  for  every  unit  charge  that 
strikes  the  wall.  It  can  also  be  seen  that  for  any  given  inner  electrode 
radius  the  intensity  was  roughly  proportional  to  the  driving  current  and 
independent  of  pressure  or  nature  of  the  gas  and  the  polarity  of  the  center' 
electrode.  There  does  appear  to  be  a  tendency  for  the  intensity  to  decrease 
as  the  electrode  radius  increases,  but  it  is  relatively  weak  and  could  easily 
be  associated  with  a  systematic  error  in  the  integration  of  the  detector  signals. 

It  is  interesting  to  note  that  the  photon  intensity  indicated  is  the  right 
order  of  magnitude  to  provide  sufficient  photoelectric  emission  from  the 
cathode  to  carry  the  primary  current.  It  is  also  sufficient  to  produce  sub¬ 
stantial  ionization  of  the  gas  ahead  of  the  current  sheet  and  we  may  reason¬ 
ably  expect  this  gas  to  have  at  least  the  room  temperature  conductivity  of 
a  fully  ionized  gas. 

Performance  without  Center  Electrode 

The  effect  of  omitting  the  conducting  center  electrode  on  the  current, 
voltage  and  speed  characteristics  of  the  shock  tube  is  shown  in  Fig.  9.  The 
records  on  the  left  were  obtained  using  a  brass  button  electrode  2.  2  cm  in 
diameter  and  1  cm  long  connected  to  a  l/4-inch  copper  rod  running  the  re¬ 
maining  length  of  the  shock  tube;  the  records  on  the  right  were  obtained 
using  the  same  button  electrode  with  the  copper  rod  replaced  by  a  30  mil 
insulated  wire  for  the  purpose  of  making  a  voltage  connection.  By  compar¬ 
ing  the  records  of  the  driving  voltage  ^ -^2  current  I^,  it  can  be  seen  that 
the  power  input  to  the  shock  tube  was  nearly  the  same  with  or  without  a  cen¬ 
ter  electrode  although  the  impedance  was  somewhat  higher  in  the  former 
case.  However,  the  fraction  of  the  power  input  which  goes  into  magnetic 
acceleration  of  the  gas  is  very  different  in  the  two  cases.  When  the  center 
electrode  is  present,  the  voltage  ^13  =  ■  ^32  ^-ssociated  with  expansion 

of  the  current  loop  is  over  97%  of  during  the  entire  test  time;  when  the 
center  electrode  is  absent  falls  rapidly  to  less  than  25%  of  'while 
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y 2,2  ^  complementary  manner.  This  indicates  that  when  the  center 

electrode  is  omitted  the  shock  tube  becomes  primarily  a  heater  and  there 
is  little  expansion  of  the  current  loop. 

A  substantial  decrease  in  the  speed  of  the  current  sheet  also  occurs 
when  the  center  electrode  is  omitted.  This  can  be  seen  from  the  UV  speed 
records  in  Fig.  9  and  the  plots  in  Fig.  10  of  the  sheet  speed  u^  and  the  flux 
speed  Uj^  as  functions  of  time.  For  purposes  of  this  comparison,  the  effec¬ 
tive  diameter  of  the  axial  current  used  in  Eq.  (2)  to  compute  Uj^  when  the 
center  electrode  was  absent  was  taken  equal  to  l/4-inch.  Since  the  depend¬ 
ence  of  Uj^  on  this  diameter  is  only  logarithmic,  the  results  would  not  be 
altered  much  by  assuming  any  other  reasonable  value  and,  in  any  case, 
the  important  qualitative  result  that  the  expansion  of  the  current  loop  is 
small  and  diminishes  rapidly  is  not  effected. 

The  operation  of  the  shock  tube  was  also  studied  under  conditions 
where  all  but  the  first  few  inches  of  the  center  electrode  was  insulated  or 
replaced  by  teflon  or  pyrex  rods,  and  was  found  to  be  substantially  the  same 
as  when  the  center  electrode  was  omitted. 

Although  the  operating  characteristics  of  shock  tubes  without  center 

electrodes  are  of  somewhat  peripheral  interest  as  far  as  the  present  paper 

is  concerned,  we  have  included  this  discussion  to  show  the  importance  of  a 

conducting  center  electrode  in  obtaining  magnetic  drive  and  because  of  the 

17-21 

similarity  between  this  geometry  and  that  of  the  "conical"  shock  tubes. 

DISCUSSION  OF  RESULTS 

It  appears  from  the  results  reported  above  that,  under  a  wide  variety 
of  operating  conditions,  only  a  small  fraction  of  the  total  drag  in  the  magnetic 
annular  shock  tube  is  due  to  the  inertia  of  the  test  gas.  In  this  section  we 
shall  explore  other  possible  sources  of  drag  to  see  if  we  can  obtain  even  a 
qualitative  understanding  of  the  results.  The  sources  considered  are:  (1) 
drag  due  to  the  flow  of  ion  current,  (2)  drag  associated  with  boundary  layer 
growth  on  the  electrodes  and  (3)  drag  due  to  the  ablation  of  material  from 
the  insulator  at  the  driver  end  of  the  shock  tube. 
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Ion  Current  Drag 

22 

It  has  been  pointed  out  by  Thom,  Norwood  and  Jaluka  that  if  the 
ions  carry  an  appreciable  fraction  f  of  the  current  in  a  magnetic- accelerator , 
there  will  be  a  drag  associated  -with  these,  ions  striking  the  cathode.  To 
calculate  this  drag  we  note  that  each  ion  car rie  s  with  it  an  axial  component 
of  momentum  mUg  and  that  the  .number  of  ions  striking  the  cathode  per  sec¬ 
ond  is  fl^/e.  The  product  of  these  factors  is  the  total  drag  due  ion  current 
and  can  be  equated  to  the  total  magnetic  pressure  to  obtain  an  expression 

for  u  .  The  result  is 
s 


Although  it  is  obvious  that  Eq.  (lo)  does  not  predict  a  limiting  speed  and 
thus  cannot  provide  the  explanation  we  are  seeking,  it  is  interesting  to  com¬ 
pare  the  magnitude  of  the  predicted  speeds  with  those  actually  observed. 
Taking  as  an  example  the  conditions  corresponding  to  the  data  of  Figs.  5 
and  6  and  assuming  f  =  1  we  find  values  for  u^  of  130,  65,  9.  3  and  6,  5  cm/jUs 
respectively,  for  H^,  He,  Air  and  Ar  as  compared  to  an  observed 
speed  of  6  cm//is.  Thus,  if  the  ions  carried  all  the  current  the  drag  could 
be  significant  for  Air  or  Ar  but  is  completely  negligible  for  or  He. 

Viscous  Drag 

We  consider  next  the  possibility  of  drag  due  to  the  growth  of  a  vis¬ 
cous  boundary  layer  on  the  electrodes  behind  the  current  sheet.  Unfortu¬ 
nately,  there  has  been  no  theoretical  treatment  of  such  a  highly  cooled 
magnetohydrodynamic  boundary  layer,  and  it  is  extremely  difficult  to 
estimate  its  thickness  due  partly  to  the  uncertainty  in  conditions  behind 

the  current  sheet  and  partly  to  the  fact  that  both  the  kinetic  mean  free 
15  3 

path  X  (cm)  =  10  /n  (cm  )  and  the  characteristic  ion  Larmor  radius 
(cm)  =  A  u^  (cm/jtls)  r  (cm)/2I  (ka)  in  the  range  of  interest  are  the 
same  order  of  magnitude  as  the  annulus  spacing.  Therefore,  our  discus¬ 
sion  of  the  viscous  boundary  layer  will  necessarily  be  limited  to  considera¬ 
tions  of  a  qualitative  nature. 

We  observe  first  that,  if  viscous  effects  are  to  provide  an  explana¬ 
tion  of  the  limiting  speed,  the  boundary  layer  behind  the  current  sheet  must 
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grow  sufficiently  to  close  the  annulus.  Otherwise  there  can  be  no  appre¬ 
ciable  flow  of  current  across  the  annulus  behind  the  leading  sheet,  and  the 
magnetic  pressure  will  be  the  same  behind  the  sheet  as  it  is  at  the  driver 
end.  For  this  to  occur  within  the  available  test  time,  the  boundary  layer 
would  have  to  grow  at  a  speed  of  at  least  ten  percent  of  the  sheet  speed, 
which  seems  extremely  unlikely  in  view  of  the  strong  cooling  that  occurs 
at  the  walls.  On  the  contrary  it  is  more  likely  that  particles  striking  the 
walls  will  tend  to  stick  and  that  the  region  behind  the  current  sheet  will  be 
cryogenically  pumped. 

If  the  boundary  layer  did  grow  sufficiently  to  close  the  annulus  during 
the  test  time,  however,  we  should  expect  an  attenuation  of  both  the  flux  and 
sheet  speeds  since  the  drag  would  increase  with  the  length  of  the  closed 
region.  As  can  be  seen  in  Figs.  5  and  6  no  such  attenuation  was  observed 
and  both  the  flux  and  sheet  speeds  were  remarkably  uniform.  We  may 
further  note  that  if  viscous  effects  were  important  they  should  be  most 
apparent  in  the  case  of  small  annulus  spacing.  This  is  a  contradiction  of 
the  data  in  Figs.  3  and  7  which  shows  that  the  observed  speeds  were  closer 
to  the  calculated  speeds  for  small  annulus  spacing  than  for  large. 

Although  a  quantitative  treatment  of  the  boundary  layer  problem  is 
to  be  desired,  we  believe  the  qualitative  arguments  given  above  are  suffi¬ 
cient  to  rule  out  any  serious  possibility  that  the  limiting  speed  observed 
can  be  the  result  of  pure  viscous  drag. 

Drag  Due  to  Ablation 

In  our  opinion  the  most  probable  explanation  of  the  limiting  speed  is 
that  it  is  due  to  the  inertial  drag  of  material  ablated  from  the  walls  of  the 
shock  tube  into  the  region  behind  the  current  sheet.  The  most  likely  sources 
of  such  material  are;  (1)  the  arc  spots  on  the  electrodes  at  the  edges  of  the 
current  sheet  and  (2)  the  insulating  wall  at  the  driver  end  of  the  tube. 

The  situation  encountered  in  the  first  case  is  basically  similar  to 
that  just  discussed,  except  the  boundary  layer  may  contain  considerably 
more  material  and  can  be  expected  to  be  somewhat  thicker.  The  same 
arguments  apply,  however,  and  it  appears  improbable  that  ablation  of  ma¬ 
terial  by  the  arc  spots  can  provide  an  explanation  of  the  observations. 
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The  situation  encountered,  in  the  second  case  is  more  interesting  and 
worth  serious  consideration.  We  observe  first  that  under  the  conditions 
existing  in  the  interior  of  the  shock  tube  we  must  certainly  expect  material 
to  be  ablated  from  insulating  wall  at  the  driver  end  of  the  shock  tube.  This 
material  will  quickly  become  ionized  and,  since  it  is  at  rest  in  a  strong 
electric  field,  a  current  will  flow  producing  both  joule  heating  and  accelera¬ 
tion  on  the  material.  Although  mo st  of  the  heat  produced  is  probably  con- 
vected  down  the  shock  tube  we  may  expect  that  a  fraction  of  it  will  return 
to  the  wall  by  conduction  and  produce  further  ablation.  Thus,  the  ablation 
process  becomes  self-sustaining  and  it  is  possible  to  conceive  of  a  situation 
in  which  material  is  continuously  generated  and  accelerated  down  the  shock 
tube  by  a  steady  or  quasi- steady  current  flowing  at  the  driver  end.  Under 
these  conditions  only  a  part  of  the  driver  current  supplied  to  the  shock  tube 
would  flow  in  the  current  sheet  and  both  the  sheet  and  flux  speeds  would  be 
reduced.  It  may  be  noted  that  the  situation  is  somewhat  similar  to  that  in 
a  coaxial  arc  jet  except  that  the  flow  rate  is  controlled  by  the  current  and 
the  wall  imposes  a  boundary  condition. 

An  attempt  to  illustrate  the  situation  graphically  has  been  made  in 
Fig.  11  which  shows  the  process  unfolding  on  a  space-time  diagram  together 
with  assumed  profiles  for  the  magnetic  field  B,  electric  field  E,  material 
velocity  u  and  current  density  j  at  a  fixed  time  after  the  sheet  speed  has 
reached  a  steady  value.  The  flow  has  been  divided  into  five  regions:  (1) 
a  photo-ionized  region  ahead  of  the  current  sheet,  (2)  a  relatively 
thin  region  in  which  the  gas  swept  up  by  the  current  sheet  is  accelerated, 

(3)  a  region  containing  ablated  material  which  grows  with  time  but  in  which 
B,  E  and  u  are  constant  and  j  is  zero,  (4)  a  region  of  length  H  in  which  the 
ablated  material  is  accelerated  by  a  current  at  the  wall  and  (5)  a  thermal 
boundary  layer  of  thickness  6  from  which  heat  is  conducted  to  the  back  wall 
to  produce  ablation. 

A  full  theoretical  analysis  of  the  flow  based  on  the  model  presented 
above  is  beyond  the  scope  of  the  present  paper.  However,  we  can  obtain 
some  insight  into  the  model  by  treating  the  problem  in  reverse  and  using 
the  experimental  observations  in  conjunction  with  the  conservation  equations 
to  obtain  an  estimate  of  the  conditions  which  must  exist  in  the  thermal  layer. 
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Our  analysis  will  be  limited  to  order  of  magnitude  considerations  and  we 
shall  neglect  all  radial  gradients.  .  We'  shall  also  assume  that  the  thickness 
of  the  thermal  layer  is  small  compared  to  the  distance  in  which  the  flow  is 
accelerated,  and  that  the  flow  is  at  least  quasi-steady  so  the  partial  deriva- ' 
tives  with  respect  to  time  may  be  neglected.  Under  these  conditions  the 
equations  for  conservation  of  mass,  momentum  and  energy  for  the  thermal 
layer  are 

m  ntu-  -  m  (11) 

a  o  0  a 

Ap  +  m  ■  u-  =  j  B  6  (12). 

^  a  0  o  o 

m  (u^  +  U;.^/2)  =  j  E  6  (13) 

a  '  a  .  0  '  “^o  o  .  '  ' 

where  m  is  the  mean  molecular  weight  of  the  ablated  material,  m  is  the 

Si-  ^ 

mass  ablation  rate,  Ap  is  the  pressure  rise  across  the  thermal  layer,  u 

ci 

is  the  energy  required  to  ablate  a  unit  mass  of  material  and  the  subscripts 
o  and  6  denote  quantities  evaluated  at  the  wall  and  the  edge  of  the  thermal 
layer,  respectively.  In  obtaining  Eq.  (13)  we  have  set  the  term  (uh-(cX)  dh/dz)g 
representing  the  net  enthalpy  transport  by  convection  and  conduction  across  the 
boundary  between  regions  4  and  5  in  Fig.  11  equal  to  zero.  This  defines  the 
thickness  of  the  thermal  layer 


6  =  Xg/Mg  (14) 

where  M  =  u/c  is  the  Mach  number  of  the  flow  and  X  is  the  mean-free-path 
for  energy  transport. 

If  the  density  ratio  across  the  current  sheet  is  large  then  the  thick¬ 
ness  of  the  slug  of  test  gas  will  be  small  and  the  final  velocity  of  the  abla¬ 
tion  products  -will  be  approximately  equal  to  the  speed  of  the  current  sheet. 
Under  these  conditions  the  characteristic  length  £  in  which  the  ablation 
products  are  accelerated  may  be  defined  by  the  equation 


B  /  =  m  u 


'o  o  as 

and  the  overall  equation  for  momentum  balance  is 


^  z 

B  ^/Z  U  =  n,m,u  +  m  u 
o'  11s  as 


(15) 

(16) 
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.  where  rij^  and  rtij^  are,  ■  respectively,  the  number  density  and  molecular 

weight  of  the  gas  ahead  of  the  leading  current  sheet.- 

In  addition  to  the  conservation  equations,  we  require  an  Ohms  law.' 

We  anticipate  that  the  gas  in  the  thermal  layer- will  be  only  slightly  ionized 

and  that  ion  slip  will  play  an  important  role  in  determining  the  conductivity. 

-  23 

■The  appropriate  Ohms  law  is 

j= - -  (E.+  uB).  (17) 

(1  +  0)  T  CO.T  .) 
e  e  1  1 

2 

where  a  =  n^e  T^/m^  is  the  gas  conductivity,  U)  is  the  cyclotron  frequency, 
T  is  the  collision  time  and  the  subscripts  e  and  i  refer  to  electrons  and  ions, 
respectively. 

To  obtain  an  explicit  solution  to  our  problem  Eqs.  (11)  through  (17) 
must  be  combined  with  a  set  of  rate  equations  describing  the  various  chemi¬ 
cal  and  physical  processes  which  together  determine  the  internal  state  of 
the  ablation  products.  Since  we  are  not  currently  in  a  position  to  even  con¬ 
sider  doing  this,  we  propose  instead  to  use  these  equations  to  deduce  the 
conditions  which  must  exist  in  the  thermal  layer  from  the  experimental 
observations.  To  avoid  algebraic  complications  we  shall  assume,  subject 
to  later  verification,  that 


A  «m  u-  , 
p  a  0 

(18a) 

u^./2  «u^, 

0'  a 

(18b) 

l«a?  T  OJ.T. 

(18c) 

e  e  1  1 

«E  /B 

0  o  o 

(18d) 

Under  these  conditions,  Eq.  (11)  to  (17)  may  be  combined  to  give 
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1 

Ml 


(^)(^) 


(20) 


(21) 


(t)=(^)(^)(t) 


(22) 


where  Uj^  is  defined  by  Eq.  (4),  E^/B^  has  been  set  equal  to  and  Q  and 

are  respectively  the  cross-sections  for  energy  and  momentum  transport. 

We  may  now  use  the  data  of  Figs.  3  and  7  to  evaluate  Eqs .  (19) 

through  (22).  We  assume  tha,t  for  lucite  which  was  the  material  of  our 

insulator  u  ~  1  cm/fXs  corresponding  to  an  effective  heat  of  ablation  of 
7  ^ 

4x10  joules/kg  and  that  the  heavy  particle  temperature  of  the  dissociating 
ablation  products  in  the  thermal  layer  will  be  in  the  vicinity  of  10,  000°K 
corresponding  to  c  ~0.  5  cm/fis  for  particles  having  a  mean  weight  of 
10  atomic  units.  Then  for  Uj^  ~  20  cm/jis,  Uj^  ~  3  cm/jds  and  Ug~  6  cm/jLl s 
which  are  typical  values  for  the  velocity  limited  region,  we  find  from 
Eqs.  (17),  (19)  and  (22) 


m  /n,  m-  u  1 0 

a'  1  1  s 

(23) 

6//  ' 

-  l/l8 

(24) 

(ne/n)5  - 

-  1/200 

(25) 

where  we  have  assumed  Q/Q^^  =  l/lO.  Equation  (23)  shows  that  for  the  con¬ 
ditions  chosen  the  mass  flux  of  ablation  products  is  an  order  of  magnitude 
larger  than  the  mass  flux  incident  on  the  current  sheet  while  Eqs.  (24)  and 
(25)  show  that  the  thickness  of  the  thermal  layer  relative  to  the  acceleration 
length  and  the  degree  of  ionization  at  the  edge  of  the  thermal  layer  are  both 
small  as  assumed.  We  see  further  from  Eq.  (21)  that  the  Mach  number  of 
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the  flow  at  the  edge  of  the  thermal  layer  is  close  to  unity  as  might  be  ex¬ 
pected  since  this  is  the  condition  which  determines  when  a  disturbance  can 
propagate  toward  the  back  wall. 

An  estimate  of  the  absolute  thickness  of  the  thermal  layer  for  test 

gas  densities  in  the  range  10^^  to  10^^  cm  ^  used  in  the  experiments  may 

- 1 5  2 

be  obtained  by  combining  Eqs.  (20)  and  (21).  Assuming  Q  ~10  cm  and 

—  2—3 

(mj^/m^)  ~1  we  find  6  ~10  to  10  cm  which  is  very  small  compared  to 
typical  shock  tube  dimensions.. 

We  may  now  check  the  validity  of  assumptions  (18a)  to  (18d)  made 
in  obtaining  these  results.  It  can  be  seen  immediately  from  Eq.  (19)  that 
(18b)  and  (18d)  are  both  satisfied  for  the  assumed  conditions.  To  determine 
whether  (17c)  is  satisfied  we  must  evaluate 

0  0 

- 1  3 

where  r  =  2.  8  x  10  cm  is  the  classical  electron  radius,  c  is  the  elec- 
o  e 

tron  thermal  speed  and  Q  is  the  Coulomb  cross-section.  Inserting  the 

®  4  5/5 

numerical  values  used  above  in  Eq.  (26)  gives  CO  T  CO.T.  ~10  to  10  (T  /lO  ) 

6  C  1  1  6 

for  the  density  range  covered  by  the  experiments.  Thus  for  T^  »  300  K, 
CO^T^CO^T^  »  1  and  assumption  (17c)  is  satisfied. 

The  final  condition  which  must  be  met  to  insure  the  consistency  of 
our  analysis  is  (17a).  To  determine  when  this  is  the  case,  we  note  that 
heat  conduction  to  the  wall  must  be  balanced  by  ablation  so  that 

rh  u  ^  ~  (c  X)  / -^2.  (27) 

a  a  6  6 

where  y  is  an  effective  ratio  of  specific  heats  for  the  ablated  material.  It 
follows  from  Eqs.  (19),  (21)  and  (27)  that  for  the  assumed  conditions 

which  is  of  order  unity  for  reasonable  values  of  y.  Thus  assumption  (17a) 
is  not  well  satisfied  and  the  Ap  term  in  the  momentum  equation  should  be 
retained  in  a  more  careful  analysis.  The  order  of  magnitude  of  our  results 
should  be  correct,  however. 


-19- 


CONCLUDING  REMARKS 

In  the  light  of  the  above  analysis  it  appears  likely  that  ablation  of 
the  insulator  at  the  driver  end  of  the  shock  tube  is  responsible  for  limiting 
the  speed  of  the'  current  sheet  in  the  magnetic  annular  shock  tube  employed 
in  these  experiments.  To  obtain  a  definitive  check  of  this  hypothesis  it 
would  be  desirable  to  investigate  the  effect  of  changing  the  material  of  the 
insulating  walls.  Unfortunately,  this  could  not  easily  be  done  in  conjunc¬ 
tion  with  the  present  series  of  experiments  because  the  importance  of 
ablation  was  not  fully  appreciated  until  after  the  experimental  apparatus 
had  been  dismantled.  It  would  also  be  desirable  to  conduct  a  more  thorough 
theoretical  investigation  of  the  problem  to  determine  whether  the  steady 
state  situation  postulated  can  actually  exist  and  what  factors  control  the 
ablation  rate  of  the  insulating  walls.  In  this  connection  it  may  be  noted 
that  the  agreement  between  calculated  and  measured  shock  speeds  obtained 
in  the  experiments  of  Patrick^  implies  that  conditions  can  exist  in  magnetic 
annular  shock  tubes  under  which  the  current  detaches  from  the  insulating 
wall.  Since  Patrick  uses  the  same  material  for  his  insulators  as  that  used 
in  the  present  experiment,  it  appears  that  these  conditions  must  be  related 
either  to  the  presence  of  bias  fields  or  the  methods  of  preionization  and 
initiation  of  the  discharge.  It  is  also  possible  that  the  complicated  current 

367 

and  flow  patterns  ’  ’  ‘  occurring  in  shock  tubes  having  a  large  radius  ratio 
have  an  important  effect  on  the  postulated  ablation  mechanism. 

In  concluding  we  should  like  to  point  out  that  the  general  problem  of 
a  current  flowing  parallel  to  an  insulating  wall  is  encountered  in  a  variety 
of  devices  other  than  coaxial  accelerators  including  MHD  generators,  cir¬ 
cuit  breakers  and  confined  arcs  and  a  detailed  analysis  of  the  postulated 
boundary  layer  would,  therefore,  be  of  considerable  practical  as  well  as 
theoretical  interest. 

ACKNOWLEDGEMENT 

The  author  would  like  to  express  his  appreciation  to  Sumner  Freshman 
who  was  responsible  for  operation  of  the  shock  tube  and  collected  most  of  the 
data  presented  in  this  report. 


-20- 


VOLTAGE  TAPS 
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Fig.  1  Schematic  diagram  of  magnetic  annular  shock  tube. 
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Typical  oscillograms  showing  the  current  and  voltage  waveforms 
for  positive  and  negative  operation  of  the  magnetic  annular  shock 
tube.  The  voltage  subscripts  refer  to  the  terminals  shown  in  Fig. 


Ijj/27rri  (ka/cm) 

Terminal  voltages  V32  as  a  function  of  current  per  unit  circumfe 
ence  of  the  inner  electrode  The  points  are  coded  as  in 

Fig.  3.  The  voltages  are  for  negative  operation;  voltages  for  po 
tive  operation  were  similar. 
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P,  =  5  0  /i  H  g  ,  r,  =  1 . 1  c  m 

Fig.  5  Representative  oscillograms  showing  the  output  of  the  5  tungsten 
UV  detectors  (lower  traces)  along  with  the  terminal  voltage  V32 
(upper  traces)  for  discharges  in  He,  Air  and  Ar  all  under  the 
same  conditions. 
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Fig.  6  Speed  us  of  current  sheet  as  a  function  of  time  for  conditions  cor¬ 
responding  to  those  of  Fig.  5.  Data  for  both  positive  and  negative 
center  electrode  are  shown.  Also  shown  is  the  flux  speed  uj_^  de¬ 
duced  from  the  flux  loop  voltage  using  Eq.  (Z)  and  the  parametric 
speeds  u;^,  U2  and  u  defined  by  Eqs.  (4),  (5)  and  (6). 
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Fig.  8  Intensity  of  photons  with  energies  above  the  tungsten  photoelectric 
threshold  (10  ev)  coming  from  the  arc  spot  on  the  inside  electrode 
plotted  as  a  function  of  the  number  of  unit  charges  striking  the 
electrode.  The  points  are  coded  as  in  Fig.  3.  The  intensities 
shown  are  for  a  negative  center  electrode;  the  intensities  for  a 
positive  center  electrode  were  the  same  within  the  experimental 
scatter . 
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Fig.  9  Oscillograms  showing  current,  voltage  and  speed  histories  for  a 
magnetic  annular  shock  tube  with  a  conducting  center  electrode 
(right)  and  without  a  conducting  center  electrode  (left).  The  records 
shown  are  for  a  positive  center  electrode;  operation  with  a  negative 
center  electrode  was  substantially  the  same. 
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Fig.  10  Sheet  speeds  Us  and  flux  speeds  ul  calculated  from  the  data  shown 
in  Fig.  9. 
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